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1. Introduction 

The Alpha Magnetic Spectrometer (AMS-02) onboard the International Space Station is the 
state-of-the-art cosmic ray (CR) detector. After 30 months of data acquisition, the AMS-02 Collab¬ 
oration has recently released the energy spectra in the 0.50—500 GeV range of the electron (e~), 
positron ( e + ) and inclusive ( e + + e~) fluxes and of the positron fraction {e + /{e + + e - )) [1, 2, 3]. 

For energies larger than 30 GeV, the e + e energy spectrum can be perfectly reproduced 
by a power-law dN/dE °< E a with a = 3.170 ±0.011. In the same energy range, an analogous 
power-law behaviour, but with a slope of about —2.8, characterises the positron spectrum. The 
difference in the slopes of the two spectra results in an increasing positron fraction for energies 
larger than « 10 GeV. This rise in the e + / (e ' + e~) was first indicated by HEAT [4], then measured 
by PAMELA [5] and Fermi- LAT [6] and finally confirmed with an unprecedented level of precision 
by AMS-02 [7], 

The behaviour of the positron fraction in the high energy window is inconsistent with the tra¬ 
ditional hypothesis that the positron flux is of purely secondary origin, i.e., generated by spallation 
reactions of primary cosmic rays (mainly protons and helium nuclei) impinging on the nuclei that 
populate the interstellar medium (mostly hydrogen and helium). In fact, the energy spectrum of 
secondary e + has a power-law slope of about —3.5 [8] which would give rise to a positron fraction 
that decreases with energy with a —0.3 slope, in strong tension with the observed slope which is 
about ±0.3. This excess has been interpreted in terms of an extra-component of primary e + from 
Pulsar Wind Nebulae (PWNe) (e.g. [9]), Supernova Remnants (SNRs) (e.g. [10]) or as a possible 
hint of annihilating or decaying dark matter (DM) particles (e.g. [11]). 

In Ref. [12] we have built a model of the e + emission from PWNe and SNRs and we have 
illustrated how these contributions can be invoked, together with the secondary emission, in order 
to explain AMS-02 data in their whole energy range. In Ref. [13] we, on the other hand, focus on 
the high energy region of the spectra, namely E > 10 GeV, with the purpose of inferring possible 
informations on the annihilation/decay rate of a particular class of DM candidates, that is Weakly 
Interactive Massive Particles (WIMPs). 

2. Astrophysical interpretation of AMS-02 data 

Cosmic electron and positron of astrophysical origin can be the outcome of a primary or sec¬ 
ondary process. Primary electrons and positrons are emitted by sources such as PWNe and SNRs, 
while, as already pointed out, the secondary production refers to those particles that are the product 
of spallation reactions involving primary CRs and the nuclei of the interstellar medium. In the next 
paragraphs we briefly summarize the basic features of these emission mechanisms in connection 
with the interpretation of AMS-02 data. We refer to [12, 13] for the full details of the analysis. 

2.1 Primary emission from Pulsar Wind Nebulae 

A PWN is powered by a rapidly spinning neutron star which hosts, on its surface, a strong 
magnetic field that possesses an intensity in the range of 10 7 — 10 12 G. Because of the large intensity 
of this magnetic field, charged particles can be torn away from the surface of the star. These 
particles then multiply through electromagnetic cascade effects and end up in forming the nebula, 
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that is a wind of particle and antiparticle pairs (e.g. [14]) that surrounds the pulsar. Particles of 
the nebula are accelerated up to very high energies through shock mechanisms that originate from 
the interaction of the nebula with the slow ejecta of the SNR that is associated with the pulsar. It 
is typically assumed that this acceleration process lasts for around 50 kyr. After this period, the 
nebula is disrupted and the particles are injected in the interstellar medium in a burst-like event. 

We evaluate the PWN contribution to the e ± flux by following the prescriptions outlined in 
Ref. [12]. Basically, we consider all the sources listed in the ATNF catalogue [15] and we model 
the injection spectrum of each one of them as a power-law with an exponential cut-off. Under this 
assumption, the e + flux associated to each PWN is completely described by 6 parameters: the spin 
down luminosity E, the age T, the distance d, the spectral index ypwN, the cutoff energy E c and the 
efficiency t] of the conversion of the total spin-down energy Wo into ee pairs (we address the 
reader to Ref. [12, 13] for a detailed description of these parameters). 

We assume the same cutoff energy E c = 2 TeV for every PWN, while the parameters E, d and 
T which are specific to each source, are fixed to the values reported in the ATNF catalogue. Finally, 
r] and Tpwx, which are assumed to be common to all the PWNe, are free parameters of the model 
and will be determined by fitting AMS-02 data. 

2.2 Primary emission from Supernova Remnants 

SNRs are believied to be the most powerful accelerator of CRs in the Galaxy. These sources 
are capable of accelerating electrons up to very high energies through the interaction of these par¬ 
ticles with the non-relativistic expanding shock-waves that are activated by the star explosion. As 
for the PWNe, the spectrum that is expected to be injected in the interstellar medium by SNRs can 
be described in terms of a power-law with an exponential cut-off. Informations about the e~ flux 
originated by SNRs can be inferred by assuming the radio flux B v r measured at a specific frequency 
V in the shock region to be the synchrotron emission from the electrons accelerated by the SNR. As 
widely discussed in Refs. [16, 12] this flux depends on the intensity of the magnetic field B in the 
remnant, the spectral index 7 snr< the normalization <2o of the e~ spectrum, the age T of the source 
and its distance d from the observer. 

We divide the whole population of Galactic SNRs into two categories based on their distance: 
the near SNRs are the ones for which d < 3 kpc, while the far SNRs are characterised by d > 3 
kpc. We take the parameters of the near SNRs from the Green catalogue [17]; we allow for a free 
normalization of the flux generated by the Vela SNR, which is the dominant e~ contributor among 
local sources. Clearly, this normalization reflects changes in the magnetic field Rveia- On the other 
hand, the far component is treated as an average population of sources that follow the radial profile 
derived in Ref. [18] and that share common values for <2o and 7 , which are free parameters that 
will be determined in our fit to AMS-02 data. 

2.3 Secondary production 

As discussed in Ref. [ 8 ], a crucial step in computing the secondary contribution to the e + 
fluxes, is the evaluation of the fluxes of the incoming primary CRs fficR- We address this point by 
performing a fit to the AMS-02 measurements of the proton and Helium energy spectra [19, 20]. In 
particular, we consider a fitting function given by <1>(7\' < R b reak) = Af5 p R^ Pl and <3>(R > Rbreak) = 
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Figure 1: The results of our fit to AMS-02 data of the e + +e~ (top left), e~ (top right), e (bottom left) 
and e + /(e + + e~) (bottom right) energy spectra. In each panel, the solid black line represents the best fit 
results, and it is embedded in its 3cr uncertainty band (cyan strip), while the different styles and colors of 
the other lines represent the various contributions from the component of the model: far SNRs (dot-dashed 
yellow), local SNRs (dotted green), PWNe (short dashed blue) and secondaries (long dashed red). Data from 
AMS-02 and HESS [21] are also reported. 

Pl , where R = pc/Ze is the rigidity of the nucleus of charge number Z and momen¬ 
tum p. The parameters that we obtain from the fit are: A = 26700 ±500 m 2 s 'sr^tGcV/n) -1 , 
P = 7.2 ± 0.4 and P x = 2.877 ± 0.004, P 2 = 2.748 ± 0.013 and P brea k = 220 ± 22 GV for the 
proton flux, and A = 4110 ± 80 m _2 s _1 sr _1 (GeV/n)” 1 , P = 3.5 ±0.7 and Pi = 2.793 ± 0.004, 
P 2 = 2.689 ± 0.013 and Pbreak = 187 ±18 GV for the helium flux (and for a Fisk solar modulation 
potential of 700 ± 22 MV). 

2.4 Fit results 

The electrons and positrons, that originate from the astrophysical sources discussed above, 
have to travel across the Galaxy and the heliosphere before they can reach the Earth. For the 
Galactic propagation we refer to the semi-analytical framework of the two-zone diffusion model, 
while to take into account the effects of solar modulation we adopt the force field approximation. 
The parameters related to the Galactic propagation are the ones of the MED propagation model 
[22] while the Fisk potential <p of the solar modulation is a free parameter of the fit. 

We perform a fit to the AMS-02 e~, e~ +e + , e + and e + /(e~ +e + ) energy spectra by means 
of a Markov Chain Monte-Carlo scan of the parameters space of our model performed with the 
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Figure 2: The bounds on (av) (left panel) and X (right panel) are reported for several DM annihilation/decay 
channels. 

CosmoMC package [23]. As discussed above, our model is characterised by 6 free parameters: 0, 
7 pwn- 1 ], 7 |)\vk. Qo and fiveia- The best-fit configuration that we find is associated to a reduced chi- 
squared ^“/d.o.f. = 1.03 with 191 data points and is characterised by the following values for the 
free parameters: 0 = (0.14 ±0.07) GV, ^ WN = 1.95^, r, = 3.68+™^, 7 SN r = 2.238+^^, 
Qo = (I. 23 I 0 Q 3 ) x 10 50 GeV 1 and fiveia = (30±3) pG. The best fit configuration and the 3a 
uncertainty band are displayed in Fig. 1 for e~, e +e + , e 1 and e + / (e +e + ) AMS-02 data. The 
situation illustrated in the Figure and the low chi-squared that we find are a strong indication that 
our astrophysical model provides a remarkably good explanation of AMS-02 data. 

3. Dark matter constraints 

In this Section, we expand the analysis described in Section 2.4 with the purpose of deriving 
robust constraints on the DM annihilation cross section (av) or lifetime X. We consider the 6 
annihilation/decay channels: DMDM —> {DM —?) e e . p 1 p ± x + x~, bb, it or W W . For 
each one of these channels, we determine the e ± flux by following the approach described in 
Ref. [24]: for the injected energy spectra dN/dE we refer to Ref. [25], while we asssume for the 
DM distribution in the Galaxy a NFW profile and a local density of 0.4 GeV/cnri (we have checked 
that our results and conclusions do not change significantly if other density profiles are assumed). 

The model that we use to fit AMS-02 data consists of the astrophysical contributions described 
in Sec. 2.4 with the addition of a DM component that is computed for fixed values of the DM 
mass. Therefore, the fit depends on 7 free parameters: the 6 aformentioned degrees of freedom 
that characterize the astrophysical sources plus the parameter quantifying the DM emission, that is 
(av) for the annihilating case or x if we are considering decaying DM. Within a purely bayesian 
approach, the upper limits on (av) and x are identified as the values for which the area underlying 
the posterior distributions of these parameters correspond to the desired confidence level (in our 
case, about 0.955, that corresponds to 2a bounds). 

The bounds that we obtain are displayed in Fig. 2; as it can be seen, they are particularly 
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Figure 3: Darker (lighter) shadings represent the let and 2 a contour regions that are associated to the 
best fits to AMS-02 data for a DM particle annihilating into leptonic ( left panel) and hadronic ( right panel) 
channels. In the two panels, upper limits derived from 7 -rays are shown, as discussed in the text. 

strong for the e + e~ and p ' p~ channels, for which, in the annihilating case, we are able to probe 
the thermal cross section for DM masses up to 200 and 70 GeV, respectively. 

4. Dark matter interpretation 

The purpose of this Section is to investigate the possible presence of a DM signal hidden in 
the data. To accomplish this task, we once again perform a fit to AMS-02 data within a model that 
consists of the joint contribution from both primary astrophysical sources, secondary production 
and from DM annihilation. As already widely discussed, the astrophysical emission is determined 
by 6 free parameters, while two parameters, namely (<jv) (or t) and tmdm> regulate the entity of 
DM contribution (we remark that, following a different approach with respect to the one outlined 
in the previous Section, we treat tmdm as a free parameter of the fit). 

Fig. 3 shows the la and 2a contour regions in the (moM- (crv)) plane for the annihilation 
channels under scrutiny (the contour for annihilation into tt is basically coincident with the one of 
the bb channel). Together with these best-fit configurations, the figure also reports for each annihi¬ 
lation channel the upper limit on (av) that can be derived from the study of the DM contribution to 
the Isotropic Gamma Ray Background (IGRB). In particular, solid lines represent conservative up¬ 
per limits derived in Ref. [26] by setting the 7 -ray emission from astrophysical sources ( e.g . Active 
Galactic nuclei, pulsars and Star Forming Galaxies) to the expected minimum, while dashed lines 
denote optimistic constraints derived in Ref. [27] by means of a more refined multi-component fit. 

An important remark that has to be made is that all the annihilation channels are able to provide 
an agreement with AMS-02 data that is better than the one that characterize the astrophysical 
interpretation presented in Sec. 2.4. At the same time, it is clear that only the e e and p 1 p 
channels are able at the same time to satisfy the 7 -ray bounds. In particular, the p + p~ is the most 
promising one, being characterised by a A ^ 2 ~ 49 with respect to the astrophysical fit and meaning 
that our model prefers with a significance of 6.9a the presence of DM. The best-fit parameters for 
this channel are /7?dm = 89 12 , 2 GeV and (av) = 8 1 3 x 10 26 cnvV 1 . 
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To test this possible DM hint, we consider a refined version of the astrophysical model de¬ 
scribed in Sec. 2.4 in which the efficiencies of the five most powerful PWNe of the ATNF catalogue 
are left free to vary in the fitting procedure. These five pulsars which can be found by means of 
the ranking algorithm described in Ref. [12], are Geminga, Monogem, J2043+2740, J0538+2817 
and B1742-30. If we fit the data without DM, the chi squared is 146, while the best-fit with the 
addition of DM is 122 in the case of /i * /i and 133 for x + x~ annihilation channel. It is worth 
noticing that, even within this model that allows more freedom to the astrophysical component, a 
DM annihilating into the p + p~ or r ! X channel is able to provide a better fit to AMS-02 data. 
In fact, we obtain Ax 2 ~ 24 for p + P-~ and ~ 13 for x 1 x annihilation channel associated to a 
significance of DM, with respect to our astrophysical model, of 4.5 and 3.2a respectively. The DM 
candidates has for p + p~ channel hjdm = 51^° GeV and (av ) = 2.8 ^ I x 10 26 cm 3 s 1 while 
for T + T~ channel hjdm = 140+ 2 GeV and (av) = 4 . 3 + 4 6 2 ° x 10 ~ 26 cm 3 s -1 . Therefore, the p + p~ 
DM candidate is fully consistent with y-ray upper limits shown in Fig. 3 while the x 1 x channel 
candidate is in a certain tension with the optimistic y-ray bounds. 

The contribution of the p p DM candidate to the positron fraction measured by AMS-02 is 
shown in the left panel of Fig. 4. It clearly stands out that the reason why the addition of a rel¬ 
atively light DM annihilating into leptons helps in improving the fit is that the positron flux that 
stems from this source provides a sizeable contribution at those intermediate energies at which the 
contribution from the PWNe of the ATNF catalog is relatively weak. In other words, DM helps in 
“filling the gap" where the secondary flux connects to the contribution from pulsars. 

With this being said, one can easily conclude that what has been presented until now as a possible 
hint of DM can also be seen as a mere consequence of the fact that the sources listed in the ATNF 
catalog constitute only a limited sample of the total number of PWNe in the sky. We illustrate this 
point by fitting AMS-02 data within the astrophysical model described in Sec. 2.4 but with an addi¬ 
tional PWN, whose distance, age and total emitted energy in pairs being free parameters of the 
fit. The best-fit configuration is given by a chi-square value of 146 (jy 2 /d.o.f. = 0.70) a source with 
a distance 0.59 +q [5 kpc and age 980 1 2 y|| kyr and is characterized by a Ay 2 which is comparable to 
the one given by DM. 

5. Conclusions 

In this work we have provided a quantitative characterization of the electron and positron 
data measured by the AMS-02 experiment. In particular, confirming the results already reported 
in Ref. [12], we have shown that the emission from primary and secondary astrophysical sources 
provide a perfectly viable interpretation of the experimental measurements. As a consequence, 
these leptonic observables can be used to impose strong constraints to the emission from DM 
annihilation or decay: in particular, for the leptonic channels, the bounds that we have obtained 
here, by means of a purely bayesian approach, are about as strong as the constraints that can 
be derived from y-rays investigations. Finally, we have illustrated how a contribution from DM 
annihilating into p + p~ with a mass «?dm = 50 GeV and an annihilation cross section (av) = 
2.8 x 10 26 cm 3 s _1 , which is thus fully compatible with y-ray upper limits, can improve the fit up 
to a 4.5a C.L. with respect to the purely astrophysical model. However, we have also illustrated 
how the addition of a pulsar with an age of about 1000 kyr and located at a distance of 700 pc 
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e + /(e + +e‘) POWERFUL PWN + DMDM -> mV 




Figure 4: In the left panel, the solid line reports the positron fraction as given by the jointly contribution 
of secondaries, PWNe and DM annihilating into the p p channel. The right panel shows the let (darker 
shading) and 2a (lighter shading) contour regions in the (distance, age) plane associated to the additional 
PWN which, together with the sources of the ATNF catalog, can provide the best-fit to data. 

from Earth could provide a comparably good lit of the experimental data even if a pulsar with such 
characteristics is not present in the ATNF catalog. Therefore, we conclude by saying that both 
DM and astrophysical sources are perfectly viable interpretations of AMS-02 measurements: in 
order to discriminate between these two alternatives, it will be crucial to dedicate future efforts on 
dedicated studies of the e ± emission from PWNe. 
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